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The Heisenberg antiferromagnet, which arises from the large U Hubbard model, 

V °. 

is investigated on the C%q molecule and other fullerenes. The connectivity of Cqo 

leads to an exotic classical ground state with nontrivial topology. We argue that 
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there is no phase transition in the Hubbard model as a function of U/t, and thus 
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Q\ , system undergoes a first order metamagnetic phase transition. We also consider the 
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the large U solution is relevant for the physical case of intermediate coupling. The 



S=l/2 case using perturbation theory. Experimental tests are suggested. 
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The Ceo molecule (buckminsterfullerene) has carbon atoms arranged like the vertices of 
a soccer ball. |]J We consider a neutral molecule. The active orbitals are one radial p-orbital 
for each carbon atom. When the long range Coulomb or on-site Hubbard repulsion is large 
compared with the nearest neighbor hopping t, the molecule has essentially one electron 
in each orbital. There is an antiferromagnetic Heisenberg spin-spin interaction between 
nearest neighbors, caused by superexchange. The exchange constant J = t 2 /A, where A is 
the energy difference between a state with one electron in each orbital, and a state with two 
electrons in one orbital and none in a neighboring orbital. Including on-site U and nearest 
neighbor V interactions, A = U — V. The derivation is similar to that of the t — J model 
used in connection with the superconducting cuprates. 

Estimates for real C 60 are that U is approximately 9 eV and t is 2 to 3 eV. [|2||| The real 
molecule is thus in the intermediate U regime, and not in the large U limit for which we 
calculate. The Cqq molecule is too large to numerically solve for the intermediate U ground 
state. There is evidence, however, that there is no phase transition for the Hubbard model 
in the Cqq geometry as a function of U/t. The spin correlations for intermediate U/t are 
then expected to be qualitatively similar, but smaller in magnitude, to those for large U/t. 
The evidence for a lack of a phase transition is as follows: For a finite quantum system 
at zero temperature, a phase transition as a function of the parameters occurs only if the 
quantum numbers of the ground state change. The only quantum numbers for this problem 
are spin S and angular momentum L (technically, what remains of L under the symmetry 
group of the icosahedron) . S = and L = in the limit U/t = 0, and probably also in 
the limit U/t — > oo. [|J The simplest (and we believe correct) hypothesis is that there is no 
phase transition as a function of U/t. It may be useful in this regard to consider the simple 
example of the two-site Hubbard model with two electrons. Mean field theory gives a phase 
transition with sublattice magnetization developing at finite U/t, which suggests that the 
large U limit is not continuously connected to the small U limit. The trivial exact solution, 
however, makes it clear that there is in fact no phase transition, that local moments develop 
continuously, and that the spin correlations of the large U limit develop continuously as U 



increases. 

All carbon atoms on the Cqq molecule are equivalent, but there are two slightly different 
bond lengths, 1.45A for the pentagon bonds and 1.40A for the non-pentagon bonds. |J The 
magnetic exchange constant for two neighboring sites on the same pentagon is J\. The 
constant connecting a site on one pentagon with a nearest neighbor on another pentagon is 
J 2- The Hamiltonian is 

p np 

H = Ji Yl °i ■ °k + -h Yl a j ■ <?fc, (1) 

<j,k> <j,k> 

where the first sum is over the 60 pentagon bonds, and the second over the 30 non-pentagon 
bonds. J 2 is expected to be slightly larger than J x , because the non-pentagon bonds are 
shorter. We first treat the Hamiltonian classically, so that a is a classical unit vector. This 
is the S — > oo limit. (We consider the quantum S = 1/2 case at the end of the paper.) 
The pentagons are frustrated, and cannot achieve a classical energy of —J\ per bond. The 
ground state of an isolated pentagon (5 spin system) has all spins coplanar and an energy of 
J\ cos(47r/5) = — . 80902 J\ per bond. One immediately obtains a lower bound for the energy 
of the entire ball, which is Ef, = 60 J\ cos(47r/5) — 30 J2, where we have used the fact that 
non-pentagon bonds cannot have an energy lower than — J 2 . 

It would appear that the ground state energy of Eq. (p]) cannot achieve the lower bound 
Ef,, because in simple trial states, connected pentagons interfere with each other and increase 
the energy. The classical ground state configuration was found numerically by minimizing 
the energy over the spin variables {^,0j}, i = 1,60. See Fig. (|l|). Surprisingly, the ground 
state energy is equal to the lower bound E^. The spin configuration is, however, nontrivial. 
The 5 spins in any given pentagon are coplanar, but the spins in a neighboring pentagon lie 
in a different plane. (A bond connects neighboring pentagons, but the spins on either end 
of the bond are precisely antiparallel, so that the spin planes need not be identical.) The 
ground state configuration, which has zero net moment, is the same for all positive J\ and 
J 2 . In addition to the obvious global rotational symmetry of the ground state, there is a 
discrete parity symmetry whereby each spin <jj — > —Oi. 



Since the spins in any given pentagon are coplanar, h for a pentagon can be denned to 
be the normal to the spin plane. Define a second vector rh for the pentagon, which is normal 
to the physical plane of the pentagon (a unit vector pointing away from the center of the 
ball). For a particular global spin rotation, the set of normals to the spin planes for the 
12 pentagons {hj} is a nontrivial permutation of the set of physical normals {rhj}. Let the 
global spin rotation be such that for the pentagon on the north pole, h\ = fh\. The pentagon 
on the north pole is surrounded by a first ring of 5 nearest neighbor pentagons, a second ring 
of 5 second nearest neighbors, and one pentagon on the south pole. Any of the 5 pentagons 
in the first ring have an ft equal to one of the rh for a pentagon in the second ring. Using the 
pentagon numbering convention of Fig. (P, hi = rhj^, with ( j(l), j(2), . . . , j(12) ) = (1, 7, 
9, 11, 8, 10, 6, 3, 5, 2, 4, 12). The solid angle subtended by the {hj} for 3 pentagons that are 
mutual nearest neighbors is 7 times larger than the solid angle subtended by their {rhj}. A 
topological skyrmion number is the number of times one sphere covers another sphere, much 
as a vortex number is the number of times one circle covers another circle. The skyrmion 
number for the field ft is 7, for both of the parity related ground states. [|5| Given that the 
Hamiltonian is so simple, this exotic ground state arises because of the connectivity of the 

The classical antiferromagnet on C12, C20, C70, and Csa has also been solved numerically. 
The truncated tetrahedron C12 is a smaller system with properties similar to Cqq. || It 
has 4 triangles and 4 hexagons. The classical lower bound for the truncated tetrahedron, 
El = l2Ji cos(27r/3) — 6J2 is achieved by the ground state, which has a skyrmion number 
1. 

A large number of C n fullerenes have been isolated. These compounds are closed 
(have the topology of a sphere), with 12 pentagons and a variable number of hexagons 
n h — (^/2) — 10. The average frustration decreases as the number of hexagons increases. 
We have calculated the classical ground state for the smallest molecule in this series, the 
dodecahedron C 2 o, for the most stable D 5h isomer of CV , and for the T d isomer of C 84 . fni 



The C20 molecule has not been synthesized, while C70 and Cg4 are produced in carbon arcs. 



The ground state energy for the dodecahedron is — 22.360680J, which does not reach 
the bound E^ = 30Jcos(47r/5) = —24. 270510 J. The pentagons interfere with each other, 
and prevent the system from reaching E^ . The skyrmion number of the ground state is 7. 
Assuming that all C70 bonds have the same coupling J, the ground state energy of Cm is 
—93. 346473 J, which is slightly higher than the lower bound E^ = [60 cos(47r/5) — 45] J = 
—93.541020 J. The skyrmion number of the C70 ground state is undefined. There are pen- 
tagons neighboring across the equator whose spin plane normals are precisely antiparallel, 
resulting in an undefined solid angle. Many (but not all) states that differ infinitesimally 
from the ground state have skyrmion number 7. The lowest-lying metastable configuration 
for C70 has a well defined skyrmion number 4. The T<f isomer of C§4 has the symmetry of a 
tetrahedron. In contrast to the above systems, the ground state spin configuration for C§4 
has a lower symmetry than that of the molecule. The ground state energy is —113. 892689 J, 
which does not reach the lower bound E^ = —114.541020 J. The ground state has skyrmion 
number 1, which interestingly is the same skyrmion number as the only other system inves- 
tigated with Tj symmetry, the truncated tetrahedron Cyi- 

We do not yet understand the C n problem well enough to predict the result for a general 
n isomer without doing the full calculation. The results on C 2 o, C 60 , C 70 , and C 84 are 
consistent with the hypotheses that (1) 6*60 is the unique fullerene that reaches the energy 
lower bound, and (2) the ground state has a nonzero skyrmion number when it can be 
defined. 

We now calculate the response of C 60 to a magnetic field by adding a term —hz-J2i <?i 
to the Hamiltonian. There is a first-order metamagnetic transition, which is surprising for 
an isotropic Heisenberg model. Usually magnetic anisotropy, arising from the spin-orbit 
interaction, is required for a metamagnetic transition. The symmetry of the ground state 
is different above and below the transition. Below h c , sites % and j that are diametrically 
opposite each other have identical spins, Oi = a r Above h c this symmetry is absent, but 
sites % and j that are mapped into each other by a rotation of it about the axis through the 



midpoint of one non-pentagon bond have spins related by (of, a v t1 of) = {—<Tj, —&], of)- 

In nonzero field, the spins in a given pentagon are not coplanar. The pentagon spin 
plane normal n is generalized so that a variable hij occupies each pentagon bond, with 
hi j ~ Oi x &j. The skyrmion number can be calculated for the field hij, and it is found that 
the skyrmion number is 7 both above and below the transition. 

The above results are for the lowest energy state at a given magnetic field. One can also 
obtain hysteresis loops by following metastable states while slowly changing the magnetic 
field, Fig. (§). A hysteresis loop can be complicated because of the large number of local 
minima, but it never encloses the origin. 

We have also investigated the response of C 12 , C 2 o, Cj , and C 84: to an external magnetic 
field. Of these, only C20 has a metamagnetic transition at which the magnetic moment is 
discontinuous. The other members of this group have a transition at which the moment M is 
continuous, but dM/dh is discontinuous. We do not understand this difference in behavior. 

The above M(h) calculations for S = 00 do not directly apply to Cqq, which has S = 1/2. 
Since S z is a good quantum number for the S = 1/2 system, it cannot have a magnetization 
that is linear in h for small h as shown in Fig. (2). The closest it can come to Fig. (2) is to 
follow in a staircase fashion, with a series of first order transitions at which S increases from 
to 1 to 2, etc., with a larger AS* jump at the metamagnetic transition. Since the transition 
from S — to 1 is at unobservably large fields (order of J or 1000 T), the spin susceptibility 
vanishes. For U » t, the orbital susceptibility falls as t 5 /C/ 4 from the five-membered rings. 
Since the physical U is not very much larger than t, there is an orbital contribution to the 
measured susceptibility. Elser and Haddon have estimated that the orbital contribution 



nearly vanishes for Cqq with U = 0. [|Tl|] The effect of nonzero U is hard to estimate for Cqq, 
but we note that for an isolated six-membered ring the orbital susceptibility for U/t = 4 is 
0.49 that for U — 0. The orbital and spin susceptibility should thus be very small for Cqq, 
which is consistent with measurements. [|12j 

We now discuss the S = 1/2 wavefunction. For S = 1/2, a in Eq. (|T]) is a 2 x 2 matrix. 
The simplest prescription to make an S = 1/2 trial wavefunction is to form a coherent state 
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|^o > that is a product of spinors. The spinor on each site j is quantized in the local +Zj 
direction given by the classical spin direction. The expectation < ip \H\ip > is equal to 
the classical energy. Two modifications to this prescription are required. The first is to add 
zero point spin fluctuations, which lower the energy and result in a wavefunction \ij)\ >. 
The second is to make the wavefunction the sum over coherent states representing all of the 
classical ground states (including global rotations and parity). This results in a trial state 
of total spin S = 0. 

The spin fluctuation energy is estimated in leading (second) order perturbation theory. 
This calculation gives an extremely accurate energy for the square lattice antiferromagnet. 
13] The energy shift for C 60 is 



AE = E K»c^ ) (2) 

3 C ° ~ € 3 

where each \j > has two adjacent spins flipped with respect to \ip >. For J\ — J 2 — 1, 
the ground state energy is shifted to E g = E + AE = -78.541 - 45.676 = -124.217. The 
correction to the classical energy is somewhat larger than that for the S=l/2 Heisenberg 
model on a square lattice. The reason is that in this case each site is only three-fold 
coordinated and the neighboring spins are not all exactly antialigned, so that the energy 
denominators are smaller. Another attractive variational trial state \ip2 > is the product 



of singlets on each non-pentagon bond. [14| This state has an energy E 2 = —90, which is 
considerably higher than E g . (The energy E 2 can also be reduced by adding fluctuations 
perturbatively.) We also calculated the moment reduction in second order perturbation 
theory on the coherent state. The local moment is reduced from 1 to .5590. This moment 



is smaller than that obtained for the square lattice. [JT3 

Due to quantum fluctuations, a spin 1/2 wavefunction does not have a unique skyrmion 
number. The wavefunctions \ipo > and \ipi > that we propose do, however, have unusual, 
nonvanishing spin-spin correlation functions even for widely separated spins on the molecule. 
Some of the correlation functions result from the fact that each pentagon tends to have a 
unique normal, which may be calculated from any of the five adjacent spin pairs in the pen- 



tagon, <7j x (Tj. Other nonzero correlation functions arise because the normals to different 
pentagons are related. We have calculated some correlation functions from exact diagonal- 
izations of the smaller model system, the truncated tetrahedron. These calculations indicate 



that both short and long range classical spin correlations survive quantum fluctuations. (To 
In contrast, the proposed local singlet state \ip2 > and the RVB state |2j have only short 
range antiferromagnetic correlations. 

Motivated in part by the cuprate superconductors, there has been a large effort in calcu- 
lating the properties of extra holes or electrons in a planar antiferromagnet using t — J and 



t — t' — J models. |16|] An electron doped crystal of Cqq molecules is also a superconductor, 



1 17 1 and it may be useful to do similar calculations in its more complicated spin field. 

In conclusion, we have investigated the low-energy magnetic properties of undoped C 60 , 
the related fullerenes C 2 o, C70, an d C 8 4, and the truncated tetrahedron C 12 in the strong 
interaction limit. We argue that there is no phase transition in Cqq as a function of Hubbard 
U, and thus the spin-spin correlation functions at intermediate U are expected to be similar 
but smaller than those for large U. In the classical approximation, the connectivity of 
the ball leads to an exotic magnetic ground state with nontrivial topology even for the 
simplest antiferromagnetic Heisenberg Hamiltonian. Cqq is the only fullerene investigated 
that reaches the energy lower bound, meaning that it has no frustration beyond that of an 
elementary pentagon. The spin correlations in Cqq may be measurable by inelastic neutron 
scattering, by magnetic x-ray scattering, or by 13 C NMR. 
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FIGURES 
FIG. 1. (a) The Cqq molecule is flattened into a plane by stretching the pentagon on the south 

pole. The spin directions are not changed. The spins in the center (north) pentagon are in the 

plane of the figure. A spin pointing up out of the plane of the figure is foreshortened, with an 

enlarged head. A spin pointing down into the plane has an enlarged tail. Non-pentagon bonds are 

dotted, (b) Perspective view of the spin arrangement on the ball. 

FIG. 2. Typical hysteresis loop, showing magnetic moment M as a function of applied field h, 
beginning and ending at h = 0. The couplings are J\ = 1 and 3<i = 1.2. 
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